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Abstract 


We simulate the collision of precursor icy moons analogous to Dione and Rhea as a possible origin for Saturn's 
remarkably young rings. Such an event could have been triggered a few hundred million years ago by resonant 
instabilities in a previous satellite system. Using high-resolution smoothed particle hydrodynamics simulations, we 
find that this kind of impact can produce a wide distribution of massive objects and scatter material throughout the 
system. This includes the direct placement of pure-ice ejecta onto orbits that enter Saturn's Roche limit, which 
could form or rejuvenate rings. In addition, fragments and debris of rock and ice totaling more than the mass of 
Enceladus can be placed onto highly eccentric orbits that would intersect with any precursor moons orbiting in the 
vicinity of Mimas, Enceladus, or Tethys. This could prompt further disruption and facilitate a collisional cascade to 
distribute more debris for potential ring formation, the re-formation of the present-day moons, and evolution into 
an eventual cratering population of planetocentric impactors. 


Unified Astronomy Thesaurus concepts: Saturnian satellites (1427); Impact phenomena (779); Planetary rings 


(1254); Hydrodynamical simulations (767) 


Supporting material: machine-readable table 


1. Introduction 


Whatever their origin, Saturn's rings appear to be young. 
Key observations made during the Cassini mission provided 
new measurements of the ring mass (less et al. 2019), the 
fraction of non-icy material in the rings (Zhang et al. 2017a, 
20170), and the extrinsic micrometeoroid flux at Saturn 
(Kempf et al. 2023). Together, these three factors constrain 
the ring age to be less than a few hundred megayears (Estrada 
et al. 2018; Durisen & Estrada 2023; Kempf et al. 2023). This 
is because the primarily icy rings (295% by mass, Doyle et al. 
1989; Cuzz & Estrada 1998; Zhang et al. 2017a, 2017b) are 
continuously subjected to micrometeoroid bombardment and 
are darkened over time, so an upper bound on the rings' 
exposure age can be determined by assuming that they began as 
pure ice (Cuzzi & Estrada 1998; Estrada et al. 2015). The present 
mass of the main rings is ~0.4 Mimas masses (1.5 х 1077 kg). 

Moreover, the rings appear to be losing mass at a remarkable 
rate (Hsu et al. 2018; Waite et al. 2018; O'Donoghue et al. 
2019), suggesting they are not only young but ephemeral as 
well. Durisen & Estrada (2023) showed that the process of 
direct deposition of micrometeoroids and the ballistic transport 
of their impact ejecta can account for these mass-loss rates, and 
Estrada & Durisen (2023) used these results to demonstrate that 
mass loading and ballistic transport due to micrometeoroid 
bombardment continue to drive the rings’ dynamical evolution 
once viscosity becomes too weak. These studies concluded that 
if the rings have been losing mass at more or less their current 
rate, then they may have begun with a mass of one to a few 
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times the mass of Mimas a few hundred megayears ago. They 
also find that, when ballistic transport is included, the rings' 
dynamical age and exposure age are similar—contrary to the 
view that these timescales can be vastly different, which would 
otherwise allow for the rings to look young but be dynamically 
old (Crida et al. 2019). 

The rings' geological youth and apparent transience are at 
odds with most current ideas about their origin, in which they 
are either primordial (Pollack et al. 1976; Canup 2010) or 
derived from the remnants of either a collisionally disrupted, 
Mimas-mass moon (Harris 1984; Charnoz et al. 2009; Dubinski 
2019) or a tidally split comet (Dones 1991; Dones et al. 2007; 
Hyodo et al. 2017)—catastrophes that are highly unlikely in the 
current solar system and hence presumed ancient. In the latter 
tidal disruption case, a comet or Centaur a few hundreds of 
kilometers in size is required to pass close enough to Saturn to be 
tidally disrupted, with the rings forming from some of the debris. 
The collisional disruption of a Mimas-mass moon requires a 
body tens of kilometers in size. Furthermore, the precursor moon 
must be kept close enough to Saturn long enough for the impact 
event to occur, although this could be facilitated by mean motion 
resonances with Enceladus and Dione (Dubinski 2019). How- 
ever, the flux of such bodies decreases sharply after the Late 
Heavy Bombardment (Zahnle et al. 2003), giving only a 
probability of occurrence of —10 ^ in the past few hundred 
megayears (Ip 1988; Lissauer et al. 1988; Dones 1991; Dones 
et al. 2007; Charnoz et al. 2009). 

Wisdom et al. (2022) recently explored a novel scenario for 
recent ring formation in which Saturn once had an additional 
moon between Titan and Iapetus. This Iapetus-mass satellite 
helped to maintain a secular resonance with Neptune that drove 
up Saturn's obliquity. Subsequent destablization of the moon's 
orbit by resonances with Titan could have led to a grazing 
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encounter with Saturn and a recent ring origin ~100 Myr ago, 
in addition to kicking Saturn out of the resonance. This would 
imply that Saturn is no longer in the resonance with Neptune, 
but confirming this may require more accurate measurements 
of Saturn's pole precession rate (Jacobson 2022). Detailed 
predictions are also yet to be made for the fate of the 
destabilized moon's rock and ice material, if it were sent on an 
elliptical orbit that crosses into the Roche limit (Dones 1991), 
or for its potentially disruptive effects on the rest of the satellite 
system—perhaps along similar lines to the collisional scenario 
we consider below. 

In addition to the rings, Saturn's current system of icy moons 
(interior to Titan) also presents some evidence of recent 
formation. Cuk et al. (2016, hereafter C16) demonstrated that if 
the present-day mid-sized satellite system were ancient, then 
recent tidal evolution through dynamical resonances would 
have excited the moons' orbits beyond what is seen today. In 
particular, the "smoking gun" for the moons not being ancient 
is Rhea's present inclination, which is an order of magnitude 
lower today than if it had survived crossing an evection 
resonance? at around 8.3 Rp, where Rp, = 60, 268 km is the 
equatorial radius of Saturn (at 1 bar; Lindal et al. 1985), and 
such inclinations are not readily damped (Murray & 
Dermott 1999; Chen et al. 2014). Rhea today is located outside 
evection, at 8.7 Къ, with an inclination of only 0°35, so it 
cannot have formed long ago and migrated outward. A very 
slow tidal evolution rate might also have been able to explain 
the present-day moons' unexcited orbits, but this would conflict 
with the apparent fast orbital evolution of Rhea (Lainey et al. 
2020; Jacobson 2022) and the intense tidal heating of 
Enceladus (Meyer & Wisdom 2007). The alternative is that 
Rhea formed (or re-formed) outside the resonance (C16). This 
implies a similarly recent, and thus perhaps shared, formation 
to the rings for at least some of these icy satellites. 

More recent work has shown that at least some moons of 
Saturn evolve much faster than predicted by the equilibrium tides 
assumed by C16 (Lainey et al. 2020; Jacobson 2022), possibly 
indicating resonance locking of the moons with gravity modes or 
inertial waves within Saturn (Fuller et al. 2016). The apparent 
presence of resonant tides weakens the constraints on the age of 
the system from the lack of a past Tethys-Dione 3:2 mean motion 
resonance found by C16, as it is no longer clear that these moons 
are currently on converging orbits. However, the observed fast 
migration of Rhea, if driven by resonance locking (Lainey et al. 
2020), still implies that an ancient Rhea would have crossed and 
had its orbit excited by the solar evection resonance only a few 
hundred megayears ago. 

Here we examine a new scenario based on the work of C16, 
which hypothesizes a recent dynamical instability in a 
precursor satellite system ~100 Myr ago. Any pair of moons 
with orbits and masses loosely comparable with present-day 
Dione and Rhea will be suddenly excited when the Rhea- 
analog outer moon encounters the evection resonance during 
migration. Both moons have their orbits excited to significant 
eccentricities (~0.1—0.2) and inclinations of several degrees. 
This paired excitation can naturally arise from the coupled 
migration of moons in mean motion resonance, from the inner 
moon catching up owing to the stronger tidal effects nearer 
Saturn, and from an initially separate outer moon being trapped 
just interior to evection for a long time, by a cycle of temporary 


© The evection resonance occurs when a satellite’s orbital precession period 
becomes commensurate with the planet’s orbital period around the Sun. 
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inward migration caused by the damping of evection-driven 
eccentricity (C16). Regardless of the precise evolution to this 
point, the common outcome of large eccentricities for both 
moons leads to orbit crossing and high-velocity collisions. 
Fragments and debris ejected by this impact could then 
reaccrete into the present-day moons, significantly erode other 
precursor moons leading to a collisional cascade, and perhaps 
even deliver material directly inside Saturn’s Roche limit to 
form or rejuvenate rings. 

This scenario does not rule out previous hypotheses for 
ancient satellite formation around Saturn, such as by accretion 
from primordial rings (Crida & Charnoz 2012, 2014). Such 
mechanisms would now simply apply to the formation of the 
precursor satellite system that was then recently destabilized 
and disrupted to (re-)form the moons we see today. 

Some aspects of the C16 scenario were investigated using 
smoothed particle hydrodynamics (SPH) and N-body simula- 
tions of two colliding Rhea-mass moons (Hyodo & Char- 
noz 2017), with results that suggested that the debris would 
reaccrete into a satellite without forming rings. However, this 
conclusion could be premature for a few reasons, including the 
following: only a single impact angle of 45° was considered as 
an example input for the longer-term models; the 
SPH simulations used a low resolution of 2 x 10° particles; 
the subsequent evolution assumed a simplified hard-sphere and 
dissipative coefficient-of-restitution model for the reaccretion 
of fragments, without the possibility of further fragmentation; 
and the presence of other precursor moons in the system and 
potential cascade interactions with them were neglected. 

Regarding the last issue of other satellites, the C16 scenario is 
based on a multiple-moon primordial system comparable to the 
current one. Such a precursor satellite system would likely have a 
broadly similar mass and architecture to the present system— 
perhaps with a somewhat larger total mass. The current systems of 
both Saturn and Uranus (and probably, at one time, Neptune; 
Goldreich et al. 1989) have families of mid-sized moons that 
extend inward close to their parent planets, so this expectation 
appears well justified. Furthermore, the initial distribution of 
moons cannot be random, since their radial locations must satisfy 
orbital stability criteria (e.g., Wisdom 1980; Lissauer 1995) that 
restrict the regions in which moons can orbit stably (e.g. 
Figure 1(b), Mosqueira & Estrada 2003b). We also can expect 
size sorting (increasing radius with distance), as it is a natural 
outcome of models that form mid-sized moons in a gaseous 
circumplanetary disk (Mosqueira & Estrada 2003а) or if the moons 
were “spun out" from viscously evolving, massive, primordial 
rings after planet formation (Crida & Charnoz 2012, 2014), while a 
collisional origin might allow some variation in size and 
composition (Asphaug & Reufer 2013). We therefore take the 
current Saturn system as a guide for the context and implications of 
impacts in a precursor one. 

Here we use SPH simulations with over two orders of 
magnitude higher resolution than Hyodo & Charnoz (2017) to 
model collisions based on C16 integrations. We find that such 
an impact can send significant debris throughout the system, 
including interior to Saturn's Roche limit, and can produce an 
extended distribution of massive fragments. These objects are 
likely to continue to fragment or erode and circularize with one 
another on future orbits, and they also have semimajor axes and 
eccentricities that will lead to crossing orbits, close encounters, 
and possible collisions with other precursor moons. This opens 
the possibility of a collisional cascade that leads to additional 
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Figure 1. The pre-impact orbits of p-Rhea, p-Dione, and their center of mass 
(COM) around Saturn in dark red, blue, and green, respectively, in х-у and x—z 
projections (where z = 0 is Saturn's equatorial plane), for the 3 km 5! scenario 
detailed in Section 2.4. The gray circle shows Saturn, to scale. The inset zoom- 
in panels show, as labeled, the moons and their velocities shortly before impact 
in Saturn's frame and in the COM frame. The bottom right panel shows the 
orientation for the SPH simulation, with the velocity at contact in the x- 
direction and an impact angle 0 (ignoring any tidal deformation), plus a 
possible extra rotation about the x-axis in this frame by an angle à, to change 
the orientation of the collision once placed in Saturn's frame. The initial 
separation for the simulation is set such that the time to impact is 1 hr. 


rubble piles and debris crossing within the Roche limit, which 
may be tidally disrupted or collisionally captured (if some 
material is already present) into rings. 


2. Methods 
2.1. Smoothed Particle Hydrodynamics Simulations 


SPH (Gingold & Monaghan 1977; Lucy 1977) is a 
Lagrangian method widely used to model planetary impacts 
and diverse astrophysical and other systems. Materials are 
represented by many interpolation points or "particles" that are 
evolved under hydrodynamical forces, with pressures com- 
puted via an equation of state (EOS), and gravity. SPH is 
particularly well suited for modeling asymmetric geometries 
and/or large dynamic ranges in density and distribution of 
material, such as the messy impacts considered here. However, 
standard resolutions in planetary simulations, using 107-109 
SPH particles, can fail to converge on even large-scale 
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outcomes of giant impacts (Genda et al. 2015; Hosono et al. 
2017; Kegerreis et al. 2019, 2022). 

Here we run simulations with up (о ~107° SPH particles, to 
resolve the detailed outcomes of collisions and to test for 
numerical convergence, using the open-source hydrodynamics 
and gravity code SWIFT’ (Kegerreis et al. 2019). For this study, 
we use a “vanilla” form of SPH plus the Balsara (1995) switch 
for the artificial viscosity. Future work will include more 
sophisticated modifications recently added to SWIFT that 
mitigate issues that can arise from density discontinuities in 
SPH (Ruiz-Bonilla et al. 2022; T. D. Sandnes et al. 2023, in 
preparation). As in previous works, we also neglect material 
strength, as gravitational stresses are expected to dominate 
(Asphaug & Reufer 2013), though, as smaller bodies are 
considered, this assumption will be valuable to test in future 
work. Nevertheless, we tentatively expect these standard 
simplifications to be comparatively innocuous for the high- 
speed collisions (relative to the mutual escape speed) 
considered here (Asphaug & Reufer 2013; Kegerreis et al. 
2020). 


2.2. Initial Conditions 


To distinguish the Dione-/Rhea-like moons we collide from 
their present-day analogs, we refer to the precursor satellites as 
*p-Dione" and “p-Rhea.” While a precursor system may have 
been larger in mass, we start by considering moons like those 
today, with masses of 1.10 x 107! kg and 2.31 x 10°! kg, 
respectively. Their detailed internal state cannot readily be 
predicted, so, for simplicity and in line with previous works, we 
model them both as differentiated into silicate cores and ice 
mantles, at an isothermal 100 K, with core mass fractions of 
60% and 40%, respectively. These materials are modeled with 
Tillotson (1962) basalt and water EOS (Melosh 2007), or with 
ANEOS forsterite and AQUA water (Haldemann et al. 2020; 
Stewart et al. 2020) for a small subset of comparison 
simulations. 

We generate the moons' internal profiles by integrating 
inward while maintaining hydrostatic equilibrium, and then 
we place the roughly equal-mass SPH particles to match the 
resulting density profiles using the stretched equal-area 
method." Before simulating the impact, we first run a brief 
2 hr settling simulation for each body in isolation, to allow any 
final settling to occur. The specific entropies of the particles are 
kept fixed, enforcing that the particles relax themselves 
adiabatically. The impact simulations are then evolved from 
1 Һг before contact to O9 hr after. As this is only a modest 
fraction of the orbital period around Saturn (~82hr), we 
neglect the effects of Saturn’s gravity during the 
SPH simulations, as discussed in Section 2.4. 


2.3. Impact Scenarios 


As detailed in C16, a pair of precursor moons can be 
suddenly excited by p-Rhea encountering the evection 


7 SWIFT (Schaller et al. 2018, 2023) is publicly available at www. 


swiftsim.com. 
8 The WOMA code (Ruiz-Bonilla et al. 2021) for producing spherical and 
spinning planetary profiles and initial conditions is publicly available with 
documentation and examples at github.com/srbonilla/WoMa, and the Python 
module woma can be installed directly with pip (https://pypi.org/ 
project/woma/). 

SEAGEN (Kegerreis et al. 2019) is publicly available at github.com/jkeger/ 
seagen, or as part of WOMA. 
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resonance, leading to high eccentricities, inclinations, and 
crossing orbits. We take two specific example outcomes of C16 
models, which result in collisions with relative velocities of 2 
and 3 Кт”. The orbital parameters are listed in Section 2.4. 
For comparison, the mutual escape speed is Vese ~ 580 ms |, 
so the impacts are at ~3.4vese and 5.1v,,.. These two orbit 
scenarios set the velocity at infinity, Və, for the impact initial 
conditions. Given the moons' small masses and high relative 
velocities, little additional speed (720 m s^!) is gained from 
self-gravity as the moons approach. 

The moons' orbits are rapidly and chaotically varying at this 
time, and their radii are a negligible fraction of their orbital 
semimajor axes, so we freely set a variety of impact angles to 
examine the range of possible outcomes at these speeds. For 
each of the two speed scenarios, we simulate collisions at eight 
impact angles (see Figure 1) from nearly head-on to highly 
grazing, with 3 = 5?—75? in steps of 10°. This angle positions 
the point of contact along a line from the center to the edge of 
the target moon. The impact scenario can also be rotated by an 
angle Ф (default 0°) about the x-axis in the SPH frame to 
position the point of contact anywhere on the target. While a 
Фф = 0 simulation is identical in ће SPH frame, the orientation 
in Saturn's frame and thus the resulting distribution of orbiting 
debris are affected. Therefore, we repeat the analysis of the 
debris with ¢ = 90°, 180°, and 270? for the primary suite of 
simulations. 

For each angle and speed scenario, the SPH particle mass, 
which controls the resolution, is set such that the number of 
particles per Мьһ + Mp; = 3.4 x 10?! kg is from 10? up to 107? 
in steps of 0.5 dex. In addition to the primary suite of p-Dione- 
p-Rhea impacts, we repeat these same simulations for p-Rhea— 
p-Rhea as considered by Hyodo & Charnoz (2017), although 
they used nonequal ratios between rock and ice. 

For a cleaner side test of how the results might scale with the 
precursor moons' masses, we also run a subset of impacts with 
p-Dione masses of (1, 2, 4, 8) x 10?! kg and corresponding 
p-Rhea masses of double these values. These moons are given 
а core mass fraction of >. We simulate these collisions for 
8 = 15°, 25°, and 35? for the four mass pairs. 

Finally, to further examine the numerical effects of our 
finite-particle models, each scenario at our second-highest 
resolution of 10’ particles is repeated four more times using a 
rotated orientation for the settled p-Rhea and p-Dione. 
Particularly for low-speed collisions, this can yield nonnegli- 
gible effects even at the high resolutions used here (Kegerreis 
et al. 2020, 2022)—an uncertainty that these tests will help to 
constrain. 


2.4. Transforming SPH Outputs into Saturn’s Frame 


The SPH simulations are run in the center-of-mass (COM) 
frame of the colliding moons, with the relative velocity at 
contact in the simulation box x-direction, as illustrated in 
Figure |. The simulation results must therefore be rotated and 
translated into Saturn’s frame to provide physically realistic 
pre-impact orbit conditions. The orbits of the resulting debris 
can then be computed, as detailed in Appendix B. 

The pre-impact orbital elements of p-Rhea and p-Dione are, 
for the 2km s^! collision, 


арв = 8.3 Rr, €pR = 0.12, ipR = 3°, 
app = 6.09 Rn, epp = 0.2, ір = 5° 25, 
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where a, e, and i are the semimajor axis, eccentricity, and 
inclination, respectively. The orbits intersect at the periapsis 
and ascending node of p-Rhea and the apoapsis and descending 
node of p-Dione. The relative velocity at this point of impact is 


Ver = [1.544, 0, 1.273] km 87, 
and the COM has orbital elements 
acom = 7.356 Rn, Ecom = 0.0068, icom = 0° 64. 
For the 3kms_' collision, the corresponding values are 
apr = 8.3 Ra, epr = 0.12, lg = 6°, 
app = 6.00 Rn, epp = 0.2, ipp = 1075, 


яа = [1.608, 0, 2.538] km 871, 
acom = 7.246 Rn, €com = 0.0059, icom = 1230. 


The COM of the particles in the SPH simulation is at the origin 
at rest, and the relative velocity at impact is in the x-direction. 
The particles’ positions and velocities are transformed into 
Saturn's rest frame by the following steps: 


1. Rotate about the z-axis such that v, rather than у, 15 in 
the x-direction. This rotation is by «1? and has а 
negligible effect, but we include it for completeness. 

2. Rotate about the y-axis such that the direction of the 
relative velocity changes from x to that of v,e. This 
corresponds to 3925 and 5876 for the two scenarios. 

3. Compute the position and velocity of the COM at the 
appropriate time from the input orbital elements, and then 
add them to the particles’ positions and velocities. 


As a sanity check, we confirm that these transformations are 
correct by analyzing the SPH initial snapshots to estimate the 
SPH moons' orbits, using the equations in Appendix B. These 
indeed reproduce the expected orbital elements for each moon, 
with slight variations caused by the shifted positions and 
velocities we imposed to yield the different impact angles. 
When the SPH simulation evolves further, the orbital estimates 
become less accurate, depending on how far material travels 
away from the COM. Using the final snapshots from our highly 
grazing impacts as examples, the semimajor axes remain close 
to the initial values, but by O9 hr after impact the extracted 
eccentricities can vary by ~10%. This, alongside the less 
predictable effects on the debris evolution from tidal forces 
from Saturn, results in up to tens of percent uncertainty in the 
mass of debris that crosses different orbits. However, this 
remains a smaller potential variation than that between different 
impact angles and moon masses (Section 3.2); regardless, it 
would not affect the overall conclusions of this proof-of- 
concept study, since significant masses of debris can still be 
delivered throughout the system. Nevertheless, this simplifica- 
tion should be redressed in future work. 


2.5. Identifying Debris Objects 


The surviving remnants and accumulated fragments pro- 
duced by the SPH simulations are identified using a standard 
friends-of-friends algorithm, where particles within a certain 
distance of each other—the “linking length,” /;;,,—are grouped 
together. At the relatively short time after impact of our 
analysis, most of these objects are still embedded in a debris 
field. This makes the precise choice of linking length somewhat 
arbitrary, because small changes in the linking length yield 
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Figure 2. Illustrative final snapshots in the SPH box frame from simulations with 107? particles for impact angles and speeds at infinity of 3 = 15°, v = 2.0 km s~ 
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(left) and 3 = 35°, væ = 3.0 km s^! (right). The top and bottom panels show the xz and x-y planes, respectively. Orange and blue show p-Rhea's core-rock and 
mantle-ice material, respectively, and yellow and purple the same for p-Dione. Animations are available at icc.dur.ac.uk/giant impacts/icy moons 2kms py.mp4 
and ... 3kms py.mp4, and with the same data rendered in 3D at icc.dur.ac.uk/giant impacts/icy moons 2kms.mp4 and ... 3kms.mp4. 


slightly different groups, in contrast with bodies in a vacuum 
with clear surfaces, which would be identified much the same 
for a wider range of linking lengths. For our primary resolution 
of 107? particles, we set Jinx = 15 km, which scales for other 
resolutions with the cube root of the particle mass. We test the 
sensitivity of our results to this choice by repeating the analysis 
for significantly lower and higher hink of + 20% (12 and 18 km 
for 107? particles; see Figure 3). Even smaller or even larger 
linking lengths start to no longer include the ice mantles and 
rocky cores of bodies in the same group, or connect together 
many incohesive fragments into the same group, respectively. 


3. Results and Discussion 


We first examine the general outcomes of collisions between 
precursor moons in the COM frame in which the 
SPH simulations were run in Section 3.1, before then 
considering the context of the full Saturn system in 
Section 3.2 (see Figure 1). We examine numerical and other- 
scenario comparisons for the primary-simulation results in 
Section 3.3, with a discussion of the limitations and speculative 
implications of the results in Section 3.4. 


3.1. Impact Outcomes 


For head-on to mid-angle impacts, a significant amount of 
ice and rock is ejected from both bodies, as illustrated for two 
example scenarios in Figure 2. The faster or closer to head-on 
the collision, the greater the disruption, the smaller the 
surviving remnants of the two original bodies (if any), and 
the further the debris is spread away from the impact point— 
both in and out of the orbital plane. 


Much of the ice is dispersed into diffuse debris, while the 
rock remains or more rapidly reaccretes into a large number of 
cohesive fragments. These disruptive impacts produce a 
relatively smooth mass distribution of many low-mass and 
fewer high-mass objects, as shown by Figure 3. Smaller and 
diffuse debris are not reliably resolved into discrete objects by 
the SPH simulations, so they are not included in the figure, and 
the detailed results below ~10'°kg should be interpreted with 
caution. As such, although the total mass remains the same in 
all simulations, the cumulative distributions in Figure 3 show 
only the mass of resolved fragments. The corresponding results 
for 2kms | impacts are qualitatively similar, but with а more 
rapid trend to fewer large fragments at grazing angles (see 
Appendix A). 

For more grazing collisions, the pre-impact moons are 
incrementally less disrupted. Less ice and much less rock is 
ejected, particularly for impact angles at which the moons’ 
cores do not intersect. The debris becomes dominated by the 
remnants of the two pre-impact bodies, with a smaller number 
of other fragments and diffuse ejecta (Figure 3). Although 
grazing collisions may not immediately deliver ring-forming or 
moon-disrupting material across the system, the orbits of the 
remnant moons are negligibly changed,'° and they will likely 
collide again in a future orbit (C16), at perhaps more head-on 
angles. 

Many fragments can be produced with masses above 10'8 
and even above 10!” kg (ie. order 0.1-1 Mimas masses), 
which could readily disrupt other precursor moons analogous 
to Tethys, Enceladus, and Mimas on crossing orbits (Leinhardt 


10 The post-impact semimajor axes, eccentricities, and inclinations of p-Rhea 
and p-Dione are changed by less than ~1% for 8 > 45°. 
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Figure 3. The cumulative mass distributions of objects produced by impacts 
with v4, = 3 km s after 9 hr, simulated with 107° SPH particles, colored by 
the impact angle. Additional material is ejected as diffuse debris or smaller, 
unresolved objects that bring the cumulative total mass in all cases to the same 
value beyond the low end of the horizontal axis shown. The faint lines above 
and below each primary result show the distributions for 20% smaller and 
larger linking lengths. 


& Stewart 2012), as discussed further in Section 3.2. The ice 
fractions of these larger fragments are typically a few tens of 
percent by mass, as more massive objects are better able to hold 
on to and reaccrete ice in addition to their rocky cores. Less ice 
is also ejected into distant debris by slower and/or more 
grazing impacts. For close to head-on collisions, fragments can 
be composed of up to equal mixes of p-Rhea and p-Dione 
material, while for mid-to-grazing angles, most large fragments 
are built from almost purely p-Rhea or p-Dione. 

These results are similar for the different linking lengths 
discussed in Section 2.5, and any minor differences are 
negligible compared with those between impact angles and 
speeds, as illustrated in Figure 3. The ice fractions are similarly 
insensitive to the linking length, with the larger and smaller 
values typically yielding a few percent more and less ice by 
mass, respectively. 


3.2. Trajectories in the Saturn System 


We now place the SPH results in the context of the Saturn 
system, following the procedure described in Section 2.4. 
Material is sent far throughout the system by head-on to mid- 
angle impacts, as shown in Figures 4 and 5. Large amounts of 
debris and fragments are placed onto orbits that would intersect 
with other satellites around the locations of those in the present- 
day system, including an extended population of debris that is 
directly heading deep into the Roche limit. This qualitative 
outcome is consistent across scenarios, with the angle and 
speed affecting, for example, the thickness and extent of the 
branches in Figure 4 (see Appendix A). We consider Saturn's 
Roche limit for water ice to be at 2.27 Rr, = 136,800 km, 
which corresponds to the outer edge of the A ring. This makes 
it a somewhat conservative value, as, depending on the porosity 
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Figure 4. The eccentricities and semimajor axes of the orbiting debris, for an 
example 8 = 35°, v4, = 3.0 km s^! impact, as illustrated in Figure 2 (see also 
Figure 12). Each point is one SPH particle with a mass of ~10'* kg. The dotted 
lines and corresponding particle colors indicate where an orbit's periapsis 
crosses the Roche limit or the orbits of other moons at their locations in the 
present-day system, or the apoapsis in the case of Titan. Other particles with 
semimajor axes inside and outside the impact point are colored in gray and 
light gray, respectively. The dark-red and dark-blue pairs of axis ticks indicate 
the orbital parameters of the pre-impact moons. 
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Figure 5. The inclinations of the orbiting debris, with a sign set by whether the 
longitude of ascending node, 0, is at or opposite the point of impact, for an 
example 8 = 35°, v = 3.0 km s^! impact (see also Figure 13), colored as іп 
Figure 4. The dark-red and dark-blue pairs of axis ticks indicate the orbital 
parameters of the pre-impact moons. 


and density of the debris, the effective limit could be higher 
(Tiscareno et al. 2013). 

As discussed in Section 1, a precursor system can reasonably 
be expected to have multiple mid-sized satellites in a broadly 
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Figure 6. The mass and composition of post-impact debris on orbits crossing 
the Roche limit or other moons, as a function of impact angle for va = 3.0 
km 5! scenarios. The dashed line shows the Roche-crossing mass for the same 
impact angle scenarios with orientation angle ф = 90°. The results for all 
simulations are provided in Table 1. The plotted values are exclusive, so, e.g., 
the mass shown for Enceladus's orbit does not include the material that also 
reaches Mimas or the Roche limit. 


similar architecture to the system today. The eccentric, high- 
velocity debris and fragments produced here could significantly 
erode other moons in a collisional cascade to distribute even 
more material throughout the system and into the Roche limit 
(C16), although this will have to be explored in detail with 
future modeling beyond the scope of this initial study. Head-on 
to mid-angle impacts send a few to a few tens of objects with 
masses above 10'’kg onto Enceladus-crossing orbits, for 
example, and an order of magnitude more for Tethys. 
Furthermore, a total mass of debris ranging from around that 
of Mimas to over that of Enceladus (order 1079-1020 Ке) 
reaches the orbits of the three inner moons. 

The mass that crosses each location depends on the impact 
angle and speed, as shown in Figure 6. Table 1 contains the 
results for all simulations. The distributions of debris 
trajectories are highly anisotropic, so the orientation of the 
impact point in Saturn's frame (0, Section 2.3) is also 
important and can significantly increase or decrease the 
crossing masses, depending on the impact angle. For these 
simulations, more than twice the present-day ring mass can be 
immediately sent inside the Roche limit, where it could begin 
to tidally and collisionally evolve toward forming rings (Dones 
et al. 2007), as discussed further in Section 3.4—even before 
considering further cascades and other material redistribution, 
or the wider parameter space of plausible impact scenarios 
beyond the few examples considered here. 

The Roche-entering material contains negligible rock, 
matching the present-day rings' nearly pure-ice composition 
(Zhang et al. 2017a, 2017b), while the material that could 
encounter Tethys and Enceladus can be several tens of percent 
rock by mass, depending on the scenario (Figure 6). The results 
for 2kms ! scenarios are qualitatively similar, with the 
crossing masses reduced by factors of a few, approaching an 
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order of magnitude at greater distances from the impact point, 
as detailed in Table 1. 

In addition to the range of possible impact speeds, the 
applicability of a resonant destablization like that demonstrated 
by C16 to any pair of loosely Dione- and Rhea-like analogs 
reminds us that more mass could readily be delivered to the 
Roche limit by, for example, a more massive pair of colliding 
satellites (as examined in Section 3.3), or an impact at a 
distance closer to Saturn. The latter could arise from a different 
orbital evolution of the outer two satellites, or perhaps from a 
collision between p-Dione and a p-Tethys analog—or even 
from an outer moon destabilized as in Wisdom et al. (2022). 
Correspondingly, less material would likely be delivered by a 
more distant collision or a lower-mass pair, although this is less 
likely than higher precursor masses given that the present-day 
moons must accrete from the remnants. 

A comparable ~10!°-10°° kg mass of debris is also sent 
outward by disruptive impacts to encounter Titan. Like the 
material reaching the Roche limit, this is typically also pure ice, 
although some rock can be delivered to Titan by close to head- 
on collisions. 

These initial trajectories for the debris are only the starting 
place for significant further evolution that must be examined 
and modeled in detail in future work. A common approach to 
estimating the evolution of collisional material is to assume that 
each piece of debris will evolve onto a circular orbit with the 
same angular momentum. However, disruptive collisions 
produce debris with a range of angular momenta and orbital 
energies, and detailed modeling beyond the scope of this paper 
will be required to constrain how much additional fragmenta- 
tion and spreading of material occurs. If the system were 
isolated—though that is not the expectation here in a precursor 
satellite system—then the mass of material with periapses 
within the Roche limit and the mass of material with equivalent 
circular orbits within the Roche limit could be interpreted as 
crude upper- and lower-bound estimates for the mass likely to 
evolve there. 

With these caveats in mind, Figure 7 shows the mass 
distribution of material with equivalent circular radius, 
eq = a(l — e°)cos(i)?. A significant mass of pure-ice mat- 
erial can circularize inside the Roche limit even without 
subsequent disruption, though not immediately as much mass 
as that of the present-day rings for this example collision. The 
most dispersive scenarios tested here can yield up to ~10'° kg 
of ice with deg < Rgocne, compared with ~5 х 1077 kg with 
periapses «Rgocne. Therefore, under the conservative assump- 
tion of perfectly angular-momentum-retaining collisions and 
for the specific scenarios simulated here, some further 
redistribution of material via cascade collisions or other 
processes could be required to form the rings. Such evolution 
could be supported or even made unnecessary by a more 
dispersive base impact scenario, such as with larger or faster 
colliding moons, as discussed in Section 3.3. 

We also note that a comparable but perhaps greater challenge 
could be faced by hypotheses of ring formation from the tidal 
disruption of an external body or a destabilized moon 
(Dones 1991; Hyodo et al. 2017; Wisdom et al. 2022), with 
significantly greater initial semimajor axes and eccentricities. 
Such scenarios could still lead to the creation of rings but lack 
the range of highly reduced orbital energies and angular 
momenta of the initial debris that is produced by the dissipative 
collisions considered here. 
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Figure 7. The cumulative mass of debris, if it evolved by collisional energy 
dissipation to circularize onto an equivalent orbit with the same angular 
momentum, as a function of radial distance from Saturn, for an example 
В = 25°, væ = 3.0 km 877 impact. The orange and blue lines show rock and ice 
material, respectively. The dotted line indicates the outer edge of the A ring 
at 2.27 Къ. 


Figure 7 also further demonstrates the lack of non-ice 
pollution in the initial ring-forming material, which for the 
other hypotheses could require a fortuitously ice-rich progeni- 
tor to match the observations (Dones 1991; Wisdom et al. 
2022; Estrada & Durisen 2023). Here, no rocky material has an 
equivalent circular orbit within the Roche limit or even inside 
of Enceladus's orbit in this example, and none inside of 
Mimas's orbit for any of our impact simulations. 


3.3. Convergence and Comparisons 


The total mass of debris that is placed onto orbits reaching 
the Roche limit and crossing the orbits of other moons 
numerically converges to within ~1% by 10’ particles for 
disruptive scenarios, up to a few tens of percent at other impact 
angles, as shown in Figure 8 for 3kms ' impacts, with 
complete results in Table 1. Low resolutions below ~10° 
particles yield unconverged masses that differ by many tens of 
percent for various impact angles. The masses for low-angle 
impacts converge at lower resolutions, apart from very close to 
head-on collisions, for which the mass crossing the Roche limit 
is lower and thus the resolution requirements for reliability are 
higher. This aligns with previous studies where the primary 
outcomes of high-speed, disruptive impacts converge relatively 
quickly, while impacts closer to the mutual escape speed can 
require much higher resolutions, as do extracted results of 
lower-mass features (Kegerreis et al. 2020). For 2km s! 
impacts, the mass crossing the Roche limit is lower, and thus 
convergence requires higher resolution, with most angles 
yielding differences of several percent between 10’ and 107? 
particles. In all cases, the ice fraction of the orbit-crossing 
debris converges to within a few percent by 1065-107 particles 
(Figure 8). Lower resolutions systematically underestimate the 
mass of rock delivered. 
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Figure 8. The mass of post-impact debris on orbits crossing the Roche limit 
and the composition of debris crossing Enceladus's orbit (as Roche-crossing 
material is pure ice at all resolutions), as functions of the numerical resolution, 
for v4, = 3.0 kms ! scenarios. The results for each impact angle аге 
normalized to that of the highest resolution to examine convergence. 


For the mass function of debris, higher resolutions are 
important both for converging on the distribution of large 
objects and for enabling the inspection of smaller ones. We find 
that the overall shape of the high-mass end of the mass function 
is typically well converged by 10’ particles. The distribution of 
smaller objects with masses below ^10! kg is still broadly 
similar for >10’ particles but matches less reliably. The highest 
resolution used here enables the probing of objects down to 
masses of about 1079 kg, corresponding to ~100 SPH particles, 
with potential numerical concerns for the analysis of smaller 
objects. 

We also tested repeat simulations with reoriented initial 
conditions (Section 2.3), which with infinite resolution should 
give identical outcomes but can produce nonnegligibly 
different results even at high resolutions (Kegerreis et al. 
2022). Fortunately, in these high-speed scenarios the differ- 
ences are minimal for the orbit-crossing masses (standard 
deviations <1%), and the variation in the debris mass functions 
is of a similar magnitude to that between different linking 
lengths shown in Figure 3. 

A potentially more significant limitation is the isolation of 
the SPH simulations from Saturn, as discussed in Section 2.4. 
Tidal effects would likely affect the separation and merging 
behavior of debris fragments, and the orbital estimates become 
less accurate. However, for the first several hours after the 
impact, the mass function is still rapidly evolving. By г-8-9 hr 
(or sooner in highly grazing scenarios), the results remain 
broadly consistent with time. Some fragments continue to 
merge, collide, and separate, but this yields far less variation 
than found between different impact scenarios. 

Replacing the impacting p-Dione with a larger p-Rhea 
produces similar results to the scenarios examined thus far. The 
mass functions of fragments follow the same trends shown in 
Figure 3 (and Figure 11), though as a symmetric collision there 
are two largest remnants rather than one, and at middle angles 


THE ASTROPHYSICAL JOURNAL, 955:137 (17pp), 2023 October 1 


not quite as many large, secondary fragments are produced. 
The total mass that crosses the Roche limit is within a factor of 
two of the p-Dione-p-Rhea results in most cases but is over an 
order of magnitude higher for close to head-on impacts (see 
Table 1). The mass crossing the orbits of other present-day 
moons is similarly comparable to the p-Dione-p-Rhea cases 
but is typically somewhat higher for head-on to mid-angle 
impacts and lower at more grazing angles. 

The above differences for p-Rhea-p-Rhea collisions are 
likely driven largely by the impact geometry, with some effects 
also expected from the slightly greater self-gravity of the 
moons and their post-impact remnants. At the same input 
impact angle, the centers of the two moons are more distant at 
the time of contact, so a smaller portion of the p-Rhea impactor 
intersects with the p-Rhea target than in the case of the smaller 
p-Dione impactor. This is compounded by Rhea's smaller rock 
fraction and lower density, such that grazing collisions can 
eject significant outer ice material but the rocky cores are less 
readily disrupted. This is reflected in the ice fractions of the 
orbit-crossing debris, which contain rock out to less grazing 
angles than in the p-Dione-p-Rhea case, as detailed in Table 1. 

We also tested a few scenarios with systematically varying 
masses of the precursor moons, with the p-Rhea mass 
simplified to double that of the p-Dione. The overall outcomes 
are highly similar to the fiducial simulations at the same impact 
angle. The mass crossing the Roche limit increases with the 
colliding mass almost proportionally and roughly follows a 
power law with exponent г-0.85-0.95 depending on the impact 
angle (Table 1). This substantiates the expectation that 
significantly more mass could be distributed throughout the 
system if the precursor moons were more massive than the 
present-day analogs considered for the primary simula- 
tions here. 

In contrast to this general consistency in outcomes that we 
find across scenarios, our results show some differences from 
those of Hyodo & Charnoz (2017, hereafter HC17), who ran 
SPH impact simulations of two Rhea-mass bodies and found 
less spread of debris through the Saturn system, with a 45? 
impact delivering no material directly to the Roche limit, for 
example. They then took 10% of the particles from the 45? 
SPH simulation as inputs for five N-body integrations, with a 
different random subset of particles in each run. Further 
fragmentation was not included, but particles could merge, 
following a hard-sphere model, which they found leads most of 
the debris to reaccrete into a single body on a timescale of a 
few kiloyears. 

Our SPH simulations differ in a few respects. One numerical 
difference is the SPH resolution used, as discussed above, 
which is over two orders of magnitude higher here. However, 
even with our test simulations at a comparable resolution to 
theirs, we still find more material on orbit-crossing trajectories 
for the same impact angle. Minor variations in the initial 
conditions preclude a direct comparison of the similar 
numerical techniques: First, for the impact point and the 
COM's orbit, HC17 set a circular, equatorial orbit with a 
3 Кт 877 relative impact velocity. This places the impact point 
~13% farther away from Saturn than here. The predicted orbits 
of the debris are highly sensitive to the transformation of the 
SPH results into the Saturn system (see Section 2.4, under- 
specified in HC17 with respect to the orientation of the moons), 
and a more distant impact point reduces the inward reach of 
debris. 
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Second, НС17 collided Rhea-mass bodies with core mass 
fractions of 4096 and 6096, compared with only 4096 in both 
p-Rhea bodies here. This greater core fraction might also have 
reduced the amount of widespread ejecta. The interior 
structures (and masses) of the moons in a precursor system 
are unknown, so future work to systematically examine the 
outcomes from collisions of different initial bodies would help 
to constrain this uncertainty. Similarly, simulations that include 
Saturn and model impacts at various distances from the planet 
would reveal the importance both of the planet's gravity and of 
the location of the point of impact. 

The conclusions that we draw from the SPH results in the 
context of the wider Saturn system also differ from HC17, as 
we speculate that collisions with other inner precursor moons 
could further support the delivery of material into the Roche 
limit, even with significantly less initial spread of debris than 
we find here. However, detailed N-body simulations that 
include both fragmentation and a full precursor satellite system 
will be required to make more robust predictions. It also may 
be that a more disruptive, less grazing impact than the 45? one 
examined by НС17 will be required for successful long-term 
distribution of material, even when fragmentation is 
accounted for. 

HC17 also estimated analytically that the ~kiloyear accre- 
tion timescale of their disk, neglecting fragmentation, would be 
shorter than the timescale of viscous spreading, which C16 had 
speculated could promote the delivery of material to inside the 
Roche limit. However, the immediate production of Roche- 
bound and Roche-circularizing debris that we find here and the 
possibilities of collisional cascades with other moons could 
avoid the requirement for a ring-forming disk to spread in this 
way. Furthermore, as also noted by C16, the first set of 
satellites that reaccrete from such a disk would likely be 
numerous and unstable. A longer timescale could be required 
for these satellites to merge into the relatively stable system we 
see today, during which the likely chaotic destablization and 
collisions could further distribute debris. 

Both our primary simulations and those of HC17 used 
Tillotson (1962) EOSs (and neglected material strength), which 
do not include phase boundaries and can prompt unrealistic 
behavior for highly shocked ejecta (Stewart et al. 2020). A full 
exploration of the importance and effects of these systematic 
uncertainties is beyond the scope of this project, but we make 
an initial test of the sensitivity of our results using a repeated 
subset of the p-Rhea—p-Rhea impacts simulated with the more 
sophisticated ANEOS and AQUA EOS (Haldemann et al. 
2020; Stewart et al. 2020). The overall behavior is similar to 
the Tillotson results, although in general the ice ejecta appears 
to be less diffusely dispersed. A greater proportion of ice also 
reaccretes into objects with masses around 1077-1075 kg, 
yielding a somewhat shallower size distribution and typical 
ice fractions in fragments ranging from ~20% to 100%, 
compared with a previous г-590-4096 for these p-Rhea-p-Rhea 
impacts. Most of the total masses of orbit-crossing debris are 
within a few tens of percent of the primary results, with some 
greater differences for the slower 2 km s! collisions (Table 1). 
The variations due to the impact angle and speed continue to be 
more important than the choices of EOS, moon masses, 
compositions, and numerical resolution, though the differences 
may become more significant for future work on the detailed 
subsequent evolution of the debris. 
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3.4. Implications and Future Work 


The outputs of simulations like these provide the inputs for 
the next key stage of modeling the evolution of the initially 
eccentric and inclined fragments and debris. Combinations of 
N-body integrations and further impact simulations can then 
begin to determine the consequences for the other mid-sized 
moons and the possible growth or rejuvenation of circularized 
rings within Saturn's Roche limit. 

Returning first to a brief discussion of the initial conditions 
for this scenario, the original hypothesis of recent reaccretion of 
the inner moons by C16 was based on the assumption of 
equilibrium tides. However, as discussed in Section 1, more 
sophisticated tidal models and the possibility of resonance 
locking (Lainey et al. 2020) do not avoid the recent encounter 
that an ancient Rhea would have had with the evection 
resonance. It is unlikely that this crossing could be reconciled 
with the relatively low-e, low-i orbit of Rhea without 
significant collisional damping, possibly even reaccretion 
(see C16), so an ancient age for Rhea remains highly unlikely. 
Additionally, a resonant tidal response offers many more 
possibilities for the system to be disrupted, as the energy that 
can be put into exciting satellite orbits by resonant tides is 
orders of magnitude larger than that provided by equilibrium 
tides. In future work, we hope to revisit the potential 
destablization of the system and expand on the input models 
that could lead to this general class of recent-collision 
scenarios, whether from evection excitation or, for example, 
the ejected outer satellite proposed by Wisdom et al. (2022). 

There are also some limitations of this work that merit future 
study. While we explored a range of impact angles for two 
plausible speeds and brief tests of other moon masses, this 
scenario allows for a diverse variety of impact conditions, 
including different sizes and compositions of the colliding 
moons in addition to the angle, speed, and distance from Saturn 
of the impact. We also neglected material strength, which is 
unlikely to dominate in such high-speed impacts but could have 
uncertain effects on the details and thermodynamics of the 
ejecta—as would the use of more sophisticated EOSs than in 
the primary suite of simulations. Furthermore, the 
SPH simulations were performed in an isolated box before 
being placed around Saturn. The presence of Saturn throughout 
the hours of the simulation could promote the separation of 
fragments via tidal forces, for example, and its absence adds 
uncertainty to the values of the distributed mass. Therefore, 
while our primary conclusion that a significant amount of 
debris can be sent throughout the Saturn system to cross other 
moons and the Roche limit should not be sensitive to these 
limitations, the precise details of the debris and of the 
parameter space of successful specific scenarios may be 
expected to change as these simplifications are addressed. 

The ejected fragments and debris will continue to interact 
with each other, with the other moons and any existing rings, 
and with Saturn’s tidal field, before the system evolves toward 
its current state. However, once the present-day moons are 
mostly in place, any debris that is still in orbit—including 
secondary and further cascade ejecta—could imprint itself in 
the crater populations seen today. The craters on Saturn’s 
moons indicate a population of at least some planeto- rather 
than heliocentric impactors (Ferguson et al. 2022) and do not 
match the distributions on the other outer planets’ satellites. It 
is likely premature to consider the immediate fragment 
distributions from this study as a crater-forming population, 
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Figure 9. The cumulative number distributions of debris objects produced by 
head-on to mid-angle 3 km 5! impacts, as also characterized by Figure 3. The 
lines show fitted power laws N ос M^, with exponents b given in the legend, for 
the mass range 10'°-2 x 1075 kg. The fits do not include the many smaller 
objects that are less robustly resolved, nor the lower numbers of even larger 
objects that are less well represented by a single power law in some cases. 
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since they can be expected to evolve significantly over time, 
but they can represent an initial approach toward future 
comparisons and constraints. 

With this caveat in mind, Figure 9 shows the cumulative 
number distributions of debris objects from low to intermediate 
impact angles at an impact velocity of 3kms ! (see also 
Figure 3). The distributions can be fitted crudely with a ~1/M 
power law, with a slight suggestion of a trend to steeper slopes 
for higher impact angles, as comparatively fewer medium to 
large objects are produced by more grazing collisions. One can 
then begin to speculate on how this population of impactors 
could affect the reaccreting moons and what the impact rate on 
these bodies might be—again neglecting the fact that on these 
elliptical orbits significant erosion and fragmentation could first 
be expected. For example, a dense, planetocentric population 
may provide a sufficient “fluence” of objects to saturate 
satellite surfaces in a relatively short period of time, perhaps 
~kiloyears, so that they may appear as cratered as ancient, 
gigayear-old terrain might be from a heliocentric population 
(Lissauer et al. 1988; Zahnle et al. 2003). Furthermore, the 
range of impact speeds will also be lower for planetocentric 
orbits, so different scaling laws for impactor size to crater 
diameter may be required. The thermal state of the evolving 
moons should also be considered, in terms of the cooling time 
that might be required before craters can be preserved without 
significant relaxation (White et al. 2017). However, the low 
gravitational binding energy of the small moons means that 
large impactors may cause more erosion than heating. 

In addition to the disruption, accretion, and eventual 
surfacing implications for Saturn's inner mid-sized moons, 
approximately an Enceladus mass of debris can be placed onto 
orbits that intersect with Titan. Titan's unique atmosphere and 
surface present many unanswered questions regarding their 
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evolution and history (Nixon et al. 2018). This includes a 
possible young age for the current methane system of a few 
tens to hundreds of megayears associated with the rapid 
conversion of methane into heavier hydrocarbons. The effects 
of this high-energy delivery of a large mass of icy debris are 
beyond the focus of this study but could open up new options 
for models of Titan's evolution. Furthermore, C16 also 
proposed that resonant interactions of Titan with potentially 
long-lived outer parts of the debris disk and a transient satellite 
that accretes there could excite Titan's eccentricity enough to 
explain its significant value (70.03) today. 

Finally, we could consider this general scenario for the 
destablization and disruption of a satellite system as a potential 
evolutionary pathway for other planets’ moons and rings. 
However, among the other giant planet satellite systems in our 
own solar system, a similar evection resonance destablization 
would likely not have occurred recently, if at all. In the case of 
Jupiter, where the evection resonance lies between Io and 
Europa,! ! the resonant configuration of the inner three Galilean 
moons has long been known to be stable (Laplace et al. 1829) 
and is likely ancient, perhaps as old as the solar system 
(Goldreich & Soter 1966). At Uranus, evection lies between 
Ariel and Umbriel, neither of which is likely to have crossed 
owing to the planet's large tidal Q (Cuk et al. 2020). Lastly, 
whether such a destablization could have happened at Neptune 
is unknown, since any preexisting satellite system can be 
expected to have been disrupted owing to the capture of Triton 
(Goldreich et al. 1989; Agnor & Hamilton 2006). 


4. Conclusions 


Dynamical instabilities in a precursor system of Saturnian mid- 
sized icy moons can lead to high-velocity collisions that scatter 
fragments and debris throughout the system (Cuk et al. 2016). 
Using SPH simulations with over two orders of magnitude higher 
resolution than previous studies, we find that a range of plausible 
impacts can deliver significant mass directly inside the Roche 
limit, with a ring-like composition of pure ice. Furthermore, more 
than a Mimas mass of material—-and even more than an 
Enceladus mass in some cases—is placed onto crossing orbits 
with present-day Mimas, Enceladus, and Tethys (and Titan), 
facilitating the possibility of a collisional cascade to further 
distribute material across the system. 

Head-on to mid-angle collisions produce large numbers of 
cohesive fragments and reaccreted objects with masses above 
1075-1077 kg (—0.1-1 Mimas masses). Most of these are—at 
9hr after impact—primarily composed of rock in a more 
diffuse field of ejected ice. Grazing impacts produce little 
debris, but the colliding moons' orbits are also negligibly 
changed, so the scenario retains the opportunity of a more 
disruptive collision on a future orbit as the satellites continue to 
cross paths. Slower collisions are less effective at distributing 
material, but they can still send small amounts of debris across 
the Roche limit, and they send significantly more that would 
intersect with other precursor moons to potentially cause 
further disruption. Impacts between two equal-sized moons 
yield similar results and trends. 

Even before considering subsequent collisions and future 
modeling of the further creation and redistribution of debris, we 


!! This information for these and the other giant and ice giant planet systems 
can be obtained from the JPL Planetary Satellite Mean Elements table at ssd. 
jpl.nasa.gov. 
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find that the example impacts considered here can place over 
5 x 1077 kg of ice onto Roche-crossing orbits (and no rock) and 
1075-1077 kg of pure ice with equivalent circular orbits inside 
the Roche limit. However, that these masses are factors of a 
few above and below the mass of the present-day rings is less 
important for this proof-of-concept study than the qualitative 
placement of significant masses onto Roche- and moon- 
crossing orbits across a range of impact scenarios—to which 
the production of immediately Roche-circularizing material is 
an encouraging bonus. Furthermore, the mass on Roche- 
crossing orbits scales close to proportionally with the mass of 
the colliding moons, so the larger moons that might be 
expected in a precursor system could more readily produce 
massive rings. 

The mass of debris entering the Roche limit converges for 
moderately high numerical resolutions in most cases, while the 
mass distribution requires at least ~10’ SPH particles to 
converge reliably on the population of mid-sized to large 
fragments, let alone smaller objects. 

We conclude that the impact of two destabilized icy moons 
18 a promising scenario for the recent formation or rejuvenation 
of Saturn's rings and reaccretion of mid-sized moons. Future 
work on the long-term evolution of the orbit-crossing debris, 
combined with further and more detailed modeling of collisions 
between both icy moons and smaller fragments, will help to 
constrain the implications of this scenario for Saturn's rings, its 
moons, their craters, and other surface environments. 
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Appendix A 
Extended Results for Other Impact Scenarios 


For comparison with the head-on to mid-angle impacts in 
Figure 2, example outcomes for more grazing collisions are 
illustrated in Figure 10. Less debris is ejected, and the mostly 
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Figure 10. Illustrative final snapshots from simulations with 10"? particles for impact angles and speeds at infinity of 8 = 45°, va = 2.0 km s^! (left) and 8 = 75°, 


væ = 3.0 km s~! (right), shown as in Figure 2. 


80) 

20] — 15 
530 = 
Е —— 3% 
e —— 45 
2.25 ST z 
^ — 15 
= 2.04 
E 
8 
= 
р 1.54 
$ 
B 
3 1.04 
O 

0.91 un T TTTITIT T TTTITIT T T ITIN T TTTTTIT T TT ere T 

1016 1017 1018 1019 1020 1021 


Mass, M (kg) 


Figure 11. The cumulative mass distributions of objects produced by impacts 
with v4, — 2 km sc shown as in Figure 3. 


intact remnants of the two moons continue on little-changed 
trajectories and so are likely to recollide on a future orbit, as 
discussed in the main text. 

The mass distributions from slower impacts are shown in 
Figure 11, for comparison with the faster ones in Figure 3. The 
general outcomes and trends are similar, but the debris 
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Eccentricity 


14 16 


Semi-major axis (Ry) 
Figure 12. The eccentricities and semimajor axes of the orbiting debris, in the 


final snapshot of an example 8 = 15°, v., = 2.0 km 57! impact, as illustrated in 
Figure 2, shown as in Figure 4. 


becomes dominated by the two remnants of the initial moons 
at a less grazing impact angle. 

The resulting Saturn-frame orbits of the debris from the 
slower, more head-on scenario illustrated in the left panel of 
Figure 2 are shown in Figures 12 and 13. The overall results 
are similar to the faster, higher-angle example of Figures 4 
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Figure 13. The inclinations of the orbiting debris in the final snapshot of an 


example 8 = 15°, v = 2.0 km s impact, as illustrated in Figure 2, shown as 
in Figure 5. 


and 5. The debris is more spread out around the impact point, 
which is representative for more head-on collisions in 
general, with somewhat less material reaching the Roche 
limit, which is representative for lower speeds in general. The 
inclinations are also significantly less high, because the initial 
moons’ material has been more thoroughly mixed, and the 
debris therefore lies on orbits closer to the near-equatorial 
one of the COM. 
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Appendix B 
Orbital Elements 


After the SPH particle data are transformed into Saturn’s 
frame, as described in Section 2.4, we estimate the orbit of each 
particle or fragment from its position, r, and velocity, v. The 
vis-viva equation yields the semimajor axis, a: 


(ВП) 


where r= |r|, v = |v|, and u = С(Мъ + m), with m the particle 
or fragment's mass and G the gravitational constant. 

The eccentricity, e, and inclination, i, are then derived from 
the specific angular momentum, h =r x v: 


уха r 
Т r 
Л, 
i = cos! (5) B2 
: (B2) 
The longitude of ascending node, €), is then 
Q = cos"! h (B3) 


[n2 + h2 | 
where the sign of Q is flipped if A, < 0. 


Appendix C 
Results Tables 


The primary results from all simulations are listed in Table 1, 
arranged in groups by subset, impact angle, and speed. The 
SPH simulation data will be made available on reasonable 
request. Table 1 is available in the online version in a machine- 
readable format. 
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Table 1 


Outcomes of Icy Moon Impact Simulations 


Teodoro and Kegerreis et al. 


Mass (kg) Ice Fraction (%) 

Te. En. Mi. Ro. Ti. Te. En Mi. Ro Ti 
24x10? 74x10" 75х109 13х109 15х109 68 100 100 100 100 
1.1 x 107° 1.2 x 109 1.9 x 1089 18x10" 22x10? 52 100 100 100 100 
50x10? 36x10'% 37x10" 44х109 18x10? 62 100 100 100 100 
18x10" 2.9 x 10!¢ 1.2 x 1055 1.1 x 10! 84 x 10!5 100 100 100 100 100 
55x10? 99x10?  Á 26x10? 233x107" 1.1 x 10!° 65 100 100 100 100 
3.5 х100 3.6 x 10! 19x10? 64 x 10!5 1.7 x 109 44 100 100 100 100 
10х109 18x10? 67х105 12x10 16x 10” 69 100 100 100 100 
8.6 x 1079 1.4 x 1055 1.0 x 10P 0.0 5.1 x 10!7 100 100 100 0 100 
78x10?  Á 23x10? 85x10 11x10 5.7 x 10/8 92 100 100 100 100 
48x10? 63x10? 41x10? 63x108  17»x10!5 40 88 100 100 100 
L1x109?  Á 37x10?  Á 15x10?  Á 24x10 69x 10'8 76 100 100 100 100 
27x10? 28x10" 25х10! 0.0 3.7 x 107 100 100 100 0 100 
13x10° Á 37x10? 10x10? 76x10" 20x 10!5 72 99 100 100 100 
85x10?  Á 64x10? 33x10? 46x10" 38x10" 34 81 100 100 100 
17x10? 41x10?  Á 18x10? 12x 10!8  18x10!5 57 93 100 100 100 
35x10?  Á 26x10? 14x10" 41x105  Á 22x10" 100 100 100 100 100 
93x10? 3 Á 33x10?  Á 39x10? 21x10" 49x10" 34 100 100 100 100 
95x10? 50х109 47x10" 12x10" 92х109 34 97 100 100 100 
93x10? 36x10% 44x10? 24x10" 41х10" 33 98 100 100 100 
60x10? 86x10? 18x10" 74x10? 92х109 45 100 100 100 100 
24x10? 35x10? 38х107 46x 10!° 13x10" 100 100 100 100 100 
L0x10! 76x10" 81х109 26x10! 24x 10'° 37 100 100 100 100 
L0x10!  Á 299x105? 48x10" 48х1095  7.1x10!* 36 100 100 100 100 
L0x10!  Á 238x105? 10x10" 29х10 28x 10'° 36 100 100 100 100 
1.1 x 10?! L1x10?  s58x10' 91x10? 14х109 38 100 100 100 100 
LIx10!  78х109 12x10 44x10 26x 10% 38 100 100 100 100 
LIx10! 88х107  Á40»x105 11x10 12x 10'° 38 100 100 100 100 
1.1 x10! 64x10" 29x10! 22x10! 46x 10% 38 100 100 100 100 
1.1 x10!  92х109 53х109 11x10! 56x 10" 39 100 100 100 100 
LIx10!  44x10P^ 0.0 0.0 0.0 39 100 0 0 0 
1.1 x 10?! 56х109 42x 105^ 11x10  11x10^ 39 100 100 100 100 
1.1 x 10?! 13x10" 36x 10% 0.0 0.0 39 100 100 0 0 
243 x10? 10х109 16x10? 88х107 1.1 x 10 70 71 99 100 96 
2.6x10° 16х109  Á 69x10? 14x10? 1.2 x 107° 59 63 67 100 95 
2.4 x 1020 13x10? 25x10! 2.9 x 10!8 1.1 x 1020 66 67 96 100 96 
1.6 х 1009 22x10? 81х1095 49х1095 91x 10% 84 99 100 100 100 
23 х 1020 10х109 21x10? 71x108 79x10? 58 71 99 100 100 
37x10? 26x10 79x10? 47x10? 10x10? 40 42 98 100 100 
2.7 x 1020 1.3 x 107 3.6 x 10!° 13x10? 9.0 x 10!° 50 68 100 100 100 
58x10? 52x10 28х109 11х109 45x10! 100 100 100 100 100 
23 x 10°  Á 75x10? 3.0~x 10! L1x10? 4.4 x 10!° 48 90 100 100 100 
63x10? 16х109  Á 77x10? 57x10? 2.5 x 10!° 31 58 78 98 100 
34x10? ^ 94x10? 45x10? 18x10? 49x 10!° 41 75 100 100 100 
19x10? 10x10! 13x10! 22x10 19x10 100 100 100 100 100 
23 х 100 55x10! 28x10? 8.8 х 10!8 1.7 x 109 47 78 100 100 100 
74x10?  18х109  Á 66x10? 44x10 66х10" 32 57 88 100 100 
52x10?  Á 66x10? 33x10? 1.4 x 109 1.6 x 10? 36 72 96 100 100 
12x10? 76x10" 11x10" 44x10?  À 82x10" 100 100 100 100 100 
86x10?  35х109 19x10? 26x10 43 x 10!8 32 88 100 100 100 
93x10? Á 82x10? 56х107 14x10" 16x10" 33 89 100 100 100 
88x10?  Á 40x10?  Á 22x10? 28x10?  235»x 10/8 32 85 100 100 100 
36x10? 34x10 29x10" 17х109 36x10" 99 100 100 100 100 
99x10?  Á 24x10? 41x10?  Á 38x10" 70x10!" 33 100 100 100 100 
L0x10!  s6x107" 71x10 27x10! 29x 10'° 37 100 100 100 100 
99x10? 25х109 41x10? 33x107 47х10" 33 100 100 100 100 
98x10? ^ Á 39x10? 29х107 16х109 10x10" 33 100 100 100 100 
1.1 x 10?! 39x105 39x10 35x10! 70x 10!5 38 100 100 100 100 
LIx10!  Á 47x10  Á 11x10 89x10  30»x10P 38 100 100 100 100 
LIx10!  Á 37x10? 31x10" 30х1095  43x10!5 38 100 100 100 100 
1.1 x 10?! 12x10? 16x10 60x10? 25х109 38 100 100 100 100 
1.1 x10! 34x10" 22x10! 21х10  46x10P 39 100 100 100 100 
1.1 x10! 29х10 0.0 0.0 0.0 39 100 0 0 0 
1.1 x10! Á 29x10" 14x10!  Á 21x10?  19x10P 39 100 100 100 100 
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Table 1 
(Continued) 
Mass (kg) Ice Fraction (96) 

Te. En. Mi. Ro. Ti. Te. En. Mi. Ro. Ti. 
LIx10! Á 37x10"  Á 47x10 22x10! 24x10" 39 100 100 100 100 
17x10? 92х107 10x10" 11x10 14х10 80 100 100 100 100 
54x10?  Á 98x10?  Á 24x10? 30x107 10x10? 67 100 100 100 100 
78x10?  Á 23x10? 81x108 10x10 55x 10/8 91 100 100 100 100 
13x10? ^ Á 36x10? 99х105 73x10" 19x 10'8 71 100 100 100 100 
93x10? 3 Á 33x10? 38x10? 21x10" 50х10" 34 100 100 100 100 
L0x10!  Á 299x105? 59х107 50x10! 90x 10'° 36 100 100 100 100 
1.1 x 10?! 1.0 x108 61x 10° 94x10? 11х109 38 100 100 100 100 
LIx10! 91х109 37х10 0.0 0.0 39 100 100 0 0 
23 x10? 10x10?  Á 15x10?  11x105 11х10 70 73 100 100 98 
23x10? 10х109 Á 21x10? 69x10? 7.8 x 10!° 59 72 100 100 100 
24x10? 74x10? 30x10? 11x10? 43x 10!° 47 91 100 100 100 
22 х 100 54x10! 28x10? 86x 10!8 1.7 x 109 48 78 100 100 100 
8.7 х 10°  Á35x10? 19x10?  Á 25x10? 43x 10/8 32 89 100 100 100 
99x10?  Á 23x10? 42x10? 239x107" 71х10" 33 100 100 100 100 
LIx10! 38х105  Á 42x10" 41x10! 62x 10'° 38 100 100 100 100 
1.1 x 10?! 3.6x10!7 2.8 x 10!6 1.3 x 10P 1.3 x 10P 39 100 100 100 100 
16x10? 13x10?  Á 12x10" 78x10? 14х10 99 100 100 100 100 
52x10? 96x10?  Á 24x10? 28x10" 99x 10!8 72 100 100 100 100 
77x10? 24x10? 76x10? 97х107 52x10! 92 100 100 100 100 
12x109? Á 37x10? 94x10 71x10"  19x10!5 68 100 100 100 100 
94x10?  Á 33x10?  Á 37x10? 21x10" 51х10" 34 100 100 100 100 
1.0 x 10?! 99x10? 61x10 49х10 91х109 36 100 100 100 100 
1.1 x 10?! L1x10? 64x10 56х109 11х109 38 100 100 100 100 
LIx10! 98х109 11х10 0.0 0.0 39 100 100 0 0 
23 x10? 10x10 15x10?  15x105 10x10? 69 75 100 100 100 
23 х 1020 10х109 21x10!  Á 68x10? 78x10? 58 74 100 100 100 
24x10? ^ Á72x10?  Á 31x10?  L1x10? 43x 10!° 47 91 100 100 100 
23x10?  Á 53x10? 29x10! 8.5 х 10!8 1.7 x 109 48 79 100 100 100 
86x10?  35х109 19x10? 25x10? 42x10! 32 90 100 100 100 
99x10?  Á 23x10? 42x10? Á 39x10" 7.3 x 10!” 33 100 100 100 100 
LIx10! 43x108 52х107 39х109 61x 10'° 38 100 100 100 100 
L1x10! 42x10" 2.7 x 10!5 0.0 1.1 x 10P 39 100 100 0 100 
16x10? 15x10?  Á 20x10" 65x10 13x10 100 100 100 100 100 
48x10? 93x10! 23х 1018  Á 29x10" 9.5 х 10! 79 100 100 100 100 
76x10?  Á 24x10? 71x10? 96x10" 52х108 93 100 100 100 100 
1.2 x10 37x10 89х105 70x10" 19х10! 69 100 100 100 100 
94x10?  Á 32x10?  Á 34x10?  Á 24x10" 52х10" 34 100 100 100 100 
1.0 x 10?! 99x10? 61x10 52x10! 99 x 10!5 35 100 100 100 100 
1.1 x 10?! 12x10? 62x10 33x10? 3.3 x10" 38 100 100 100 100 
LIx10? 7.5 x10" 0.0 0.0 0.0 39 100 0 0 0 
2.3 x10?  Á 98x10? 15x10?  18x105 10x10? 69 81 100 100 100 
23x10? 99x10? Á 22x10? 67x10? 78х109 57 76 100 100 100 
25x10?  Á 68x10?  Á 32x10? 10x10? 44x 10!° 45 97 100 100 100 
23x10?  Á 53x10? 29x10? 83х108 1.6 x 10? 48 80 100 100 100 
8.7 х 1020  Á 34x10? 18x10? 24x10? 4.1 x 10/8 32 94 100 100 100 
L0x10!  Á 22x10? 41х105 42x10" 78x10" 33 100 100 100 100 
11х10!)  Á 42x10? 48x10" 39х109 60x 10'° 38 100 100 100 100 
11х10!)  51х107 9.9 x 105^ 0.0 0.0 39 100 100 0 0 
44x10? 25х10" 0.0 0.0 72x10 100 100 0 0 100 
40x10? 97x108 22x108  Á 26x10" 94x10 100 100 100 100 100 
77x10?  Á 24x10?  Á 69x10? 93x10" 53x10! 95 100 100 100 100 
L1x109? 36x10 89х105 77x10" 20х10! 72 100 100 100 100 
94x10? 30x10! 35x108  Á 27x10" 52x10! 34 100 100 100 100 
1.0 x 10?! L1x10? 74x10" 42х109 90х109 35 100 100 100 100 
1.1 x 10?! 16x10? 53x 10!6 0.0 0.0 38 100 100 0 0 
LIx10? 2. x 10!6 0.0 0.0 0.0 39 100 0 0 0 
22x10?  Á 95x10? 16x10?  Á21x105 1.0 x 10° 69 85 100 100 100 
23x10? 94x10? Á 22x10? 68x10? 7.8 x 10!° 56 83 100 100 100 
23x10? ^ Á 67x10?  Á 32x10? 11x10? 44х10! 48 99 100 100 100 
22x10? 52х109 28x10? 84х105 16x10? 49 84 100 100 100 
8.7 х 1020  Á 32x10? 18x10? 26x10? 46х108 31 100 100 100 100 
99x10? 3 Á 23x10? 45x10? 49x10" 88x 10!7 33 100 100 100 100 
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Table 1 
(Continued) 
Mass (kg) Ice Fraction (90) 

Te. En. Mi. Ro. Ti. Te. En. Mi. Ro. Ti. 
LIx10!  50х105 58х107 53х10  77x10!5 37 100 100 100 100 
LIx10! 26х10 0.0 0.0 0.0 39 100 0 0 0 
3.6x10'% . 27 x 10" 0.0 0.0 4.1x 10? 100 100 0 0 100 
29x10? 56x10? 71x10" 0.0 3.4 x10? 100 100 100 0 100 
74x10? 26x10?  Á 67x10? 87х107 52x10! 99 100 100 100 100 
1.2 x10  Á 35x10? 84x10 78x10" 19x 10!5 70 100 100 100 100 
96x10? 28x10? 35x10? 20x10" 65х10" 34 100 100 100 100 
L0x10! 89х105 61x 10!7  Á34x10/^ 31x 10'° 35 100 100 100 100 
1.1 x 10?! 1.4 x 105 0.0 0.0 0.0 38 100 0 0 0 
1.1 x 10?! 0.0 0.0 0.0 0.0 39 0 0 0 0 
23x10? 89х109 18x10?  Á 22x10 9.9 x 10!° 66 90 100 100 100 
22x10? 92x10? 23x10? 70x10  77»x10? 58 86 100 100 100 
24x10?  Á 66x10? 33x10? 10x10? 43x10? 47 100 100 100 100 
243 x10? 48x10?  Á 27x10? 83х105 15x10? 48 92 100 100 100 
89x10?  Á 31x10? 16x10? 26x10? 42x10! 31 100 100 100 100 
L0x10! ^ 21x10? 37х105 47x10" 65x10" 34 100 100 100 100 
L1x10! 40x108 54x10" 0.0 0.0 38 100 100 0 0 
11х10! 34х109 0.0 0.0 0.0 39 100 0 0 0 
22 х 100 38x10! 41x 10!8 37х107 7.5 x 10!8 81 100 100 100 100 
2.4. x10? 53x10? 1.6 x 1079 17x10 9.6 x 10!5 85 100 100 100 100 
20x10?  À 75x10? 19x10? 17x10 83x 10/8 97 100 100 100 100 
78x10? ^ Á 75x10? 10x10? 75x10" 34x10! 60 100 100 100 100 
24x10! ^ Á24x10? 26x10? 23x10" 91х10" 52 100 100 100 100 
22x10! 50x108 41x10" 51х109 20х10" 55 100 100 100 100 
23x10! 81x10" 50х109 11х109 33х109 57 100 100 100 100 
23x10! ^ Á 69x10 55x105  Á34x10^ 20x10" 57 100 100 100 100 
42x10? 217х109  Á 12x10?  Á 36x10? 97x10? 63 71 99 100 100 
41x10?  Á 23x10? 94x10  Á 34x10? 76x10? 64 79 99 100 100 
32x10? 17х109  Á 78x10? 25x10? 56х10! 70 95 100 100 100 
5.4 x 1020 13x10? 52x10? L1x10? 2.8 x 10!° 79 98 100 100 100 
18x10! ^ ÁsS4x10? 17x10?  Á 25x10  75»x10!5 51 100 100 100 100 
22x10! ^ Á 21x10? 37x10 37x10" 11x10! 54 100 100 100 100 
23x10! ^ Á 37x10?  Á 36x10" 42х109 12х10" 57 100 100 100 100 
23x10! 36х107 28x10. 25x10? 81х10 57 100 100 100 100 
15x10?  Á 25x10? 19х105 g9x 10! 51x10 100 99 100 100 100 
90x10? 94x10? 78x10" 50х109 17x10? 100 99 100 100 100 
243 x10? 14x10? 13x10" 26x10  Á 33x10" 100 99 100 100 100 
23x10? 10х109  Á 33x10? 51х105 55х10! 94 100 99 100 100 
6.9 x 1020 13 x 10  Á 27x10?  Á 231x105? 97x10! 83 100 100 99 99 
19x10! 41x10? 54x108 46x10" 14х10! 61 100 99 99 99 
23x10? 55х109 22x10? 59х105 3.8 x 10!° 80 100 100 100 100 
17x10? 69x10!  Á 26x10? 64x10 18x10? 79 100 100 100 100 
17x10? 72x10  Á 19x10? 28x 10!8 52х10! TI 99 100 100 100 
39x10?  Á 10x10? 42x10? 79x10? 68x10? 83 100 100 100 100 
2.9 х 1020 13x10? 50x10! 86x108 29x10? 86 100 100 100 100 
39x10?  13х109 32х109 233x105 77x10! 81 99 100 100 100 
67x10? 20х109 82x10 17x10? 1.3 x 107? 89 100 100 100 100 
49x10? 24x10?" 92x10?  Á 20x10? 58x10? 91 100 100 100 100 
77 х 1020  Á 24x10? 55x10! 66х108 1.5 x 10P? 86 100 100 100 100 
97x10? 38x10 16x10 40x10? 2.6 x 107° 97 100 100 100 100 
81x10? 46x10 18x10  Á 45x10? 12x10? 96 100 100 100 100 
L6x10! 41x10?  Á 10x10? 16x10? 31х10 88 100 100 100 100 


Note. Columns list the impact angle, 8; speed at infinity, va; the mass of debris on orbits crossing the Roche limit and present-day moons; and the fraction of which is 
ice. For the base suite of simulations, each entry is followed by three more lines with the SPH frame rotated about the x-axis by ¢ = 90°, 180°, and 270° (see Figure 1). 
The "type" notes other parameters that are changed from the base scenario for the following subsets of simulations, as described in Section 2.3, within which the angle 
and speed are then varied, namely, the number of particles, N (base 1077): the body colliding with p-Rhea (base p-Dione); the EOS (base Tillotson); and the mass of 
p-Dione for the separate tests of moon masses, Mpp in units of 10?! kg, with corresponding p-Rhea mass Myr = 2 Mpp. Note that the orbit-crossing masses include 
only the mass that reaches that orbit and exclude material that also reaches further ones (e.g., the Enceladus values do not include Mimas- or Roche-crossing material). 


(This table is available in machine-readable form.) 
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